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Abstract 

Relations  between  position  and  velocity  vec¬ 
tors  at  different  points  on  a  trajectory  in  a  pure 
inverse-square  field  of  force  are  derived  without 
the  use  of  geometrical  descriptors  of  the  orbit. 

An  along-track  'minimal'  transformation  variable  is 
found,  which  permits  the  direct  integration  of  the 
equation  of  motion.  The  result  is  equally  applic¬ 
able  to  elliptic,  hyperbolic,  parabolic  and  recti¬ 
linear  trajectories.  The  relationship  between  the 
transformation  variable  and  time  constitutes  an 
archetype  of  Kepler's  equation,  conventional  forms 
of  that  equation  appearing  as  special  cases.  The 
results  allow  a  further  simplification  for  recti¬ 
linear  motion,  with  the  velocity  used  as  the  along- 
track  variable.  The  'minimalist*  approach  is  also 
applied  to  the  rendezvous  problem:  Lambert's  cele¬ 
brated  theorem  reduces  to  an  obvious  observation. 
Application  of  the  theorem  to  the  rectilinear 
trajectory  allows  the  physical  interpretation  of 
parameters  introduced  by  other  authors  through  a 
purely  mathematical  analysis. 


Introduction 

Problems  involving  the  motion  of  a  body, 
idealized  as  a  point  mass,  in  an  inverse-square-law 
force  field,  with  or  without  perturbations,  are 
usually  handled  with  the  aid  of  a  general  theory  of 
orbits,  which  has  evolved  over  several  centuries 
and  includes  a  rich  treasury  of  mathematical  tools 
for  analysing,  describing  and  calculating  both  the 
instantaneous  dynamical  properties  of  the  point 
mass  and  the  geometrical  properties  of  its  path, 
the  orbit. 

However,  the  very  antiquity  of  the  dynamical 
problem,  encountered  by  astronomers  long  before 
Newton's  discovery  of  the  inverse  square  law,  has 
established  a  tradition  of  geometrical  approaches 
to  its  solution!”^  Indeed,  the  inverse  square  law 
itself  was  originally  deduced  from  the  geometrical 
properties  of  observed  planetary  orbits.  This  tra¬ 
dition  persists  to  this  day  and  the  illustrative 
nature  of  geometrical  representation  can  be  a  help 
in  the  comprehension  of  many  features  of  the  dyna¬ 
mical  behaviour.  Nevertheless,  the  preoccupation 
with  the  geometrical  properties  of  the  orbit  can 
sometimes  obscure  the  physics  of  the  motion,  divert¬ 
ing  attention  from  the  dynamical  aspects,  such  as 
energy  and  angular  momentum,  to  difficulties  asso¬ 
ciated  with  the  geometric  descriptors  of  the  orbi¬ 
tal  conic,  such  as  the  parameters  describing  the 
orientation  of  its  axes.  The  obscuration  of  the 
problem  is  particularly  noticeable  when  the  'all- 
points*  approach  via  the  geometry  of  the  orbital 
path  is  applied  to  problems  that  really  involve 
conditions  (position  and  velocity  components)  at 
just  two  points  of  the  orbit:  an  'initial'  point, 
and  the  point  reached  at  an  arbitrary  later  (or 
earlier)  instant. 


The  'minimal*  approach  of  this  paper  to  the 
two-point  problem  has  two  main  thrusts.  First,  the 
solution  of  the  equations  of  motion  is  obtained 
without  reference  to  the  geometry  of  conic  sections, 
though  the  parameters  naturally  entering  the  mathe¬ 
matical  solution  can  Teadily  be  related,  a  pos¬ 
teriori,  to  the  well-known  geometric  descriptors. 
Secondly,  instead  of  assuming  the  solution  of  a 
second-order  differential  equation  in  terms  of 
(geometrically  based)  trigonometric  or  hyperbolic 
functions,  a  universal  solution  is  obtained  by 
direct  integration,  after  a  suitable  transformation 
of  variable. 

In  section  2  a  brief  derivation  of  the  second- 
order  differential  equation  governing  the  motion  is 
followed  by  a  first  integration,  resulting  in  the 
fundamental  first-order  equation  (the  energy 
equation)  expressed  in  terms  of  the  initial  position 
and  velocity  components.  Section  3  presents  the 
'minimal*  transformation  needed  to  make  that  equa¬ 
tion  directly  integrable.  The  relation  between 
time  and  the  along-track  coordinate  so  introduced  is 
derived  in  section  4;  it  appears  as  an  archetype  of 
Kepler's  equation.  The  special  case  of  rectilinear 
motion  allows  further  simplification  and  is  treated 
in  section  5.  This  completes  the  analysis  of  the 
problem  in  which  the  conditions  at  some  'initial' 
instant  are  known  and  those  at  a  different  time  are 
to  be  found. 

In  section  6  the  results  are  applied  to  a 
modified  two-point  problem:  the  positions  at  both 
ends  of  a  trajectory  and  the  time  of  travel  are 
given,  and  the  initial  (and  final)  velocity  to 
achieve  the  rendezvous  is  calculated;  this  is 
usually  known  as  the  Lambert  Problem  and  has  an 
extensive  literature.  Of  course,  the  'minimalist' 
approach  yields  solutions  which  are  entirely  equi¬ 
valent  to  those  achieved  using  traditional  methods: 
we  start  from  the  same  vector  differential  equation 
and  set  of  boundary  conditions;  only  the  methodol¬ 
ogy  is  different. 

Finally,  in  section  7,  the  solution  of  the 
Lambert  problem  obtained  in  section  6  is  compared 
with  the  formulations  of  Lancaster  et  al  (Ref  4, 
expanded  in  Ref  5).  The  comparison  allows  a  physi¬ 
cal  interpretation  (in  terms  of  rectilinea  ‘  ’ - 
parameters)  of  the  quantities  introduced  i  '..*f  -  on 
a  purely  algebraic  basis. 


2.  Basic  dynamics 

The  motion  of  the  point  mass  is  governed  by 
the  vector  differential  equation  of  second  order: 
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Substitution  in  (6)  yields 


where  r  is  the  radius  vector  from  the  centre  of 
attraction  and  p  is  a  constant. 

If  R  and  Q  are  the  position  and  velocity 
vectors  at“an  initial  instant  (taken  as  time 
t  -  0),  then  the  motion  is  confined  to  the  plane 
through  the  centre  of  attraction  (origin)  and  con¬ 
taining  the  vectors  R  and  g  .  To  establish  the 
'minimal'  equations  of  motion  we  have  to  define  a 
reference  system  in  that  plane;  the  obvious  choice 
is  the  pair  of  orthogonal  unit  vectors  I  and  J  , 
with  I  along  R  and  the  sense  of  J  chosen  to 
make  g  .  J  non-negative  (J  is  not  unique  when 
g  .  J  ■  0) .*  Thus  we  can  write 


which  integrates  to 


2  2  2 

/dr  \  V  R  2 v 

Idl)  *  2  T  ■  constant  • 


This  is  the  energy  integral  of  the  motion  (per  unit 
mass)  multiplied  by  2.  The  quantity  W  has  the 
dimensions  of  velocity  and  is  real,  zero  or  pure 
imaginary,  depending  on  the  total  energy  of  the 
motion  (W  or  iW  is  taken  as  non-negative) .  In 
terms  of  the  initial  position  and  velocity. 


RI  (R  >  0) 


UI  +  VJ  (V  >  0).  (2) 


For  the  description  of  the  position  r  and 
velocity  q  at  an  arbitrary  time  t  ,  polar  co¬ 
ordinates  (r,0)  are  more  convenient  than  cartesian. 
Thus  we  have 

r  ■  r(I  cos  0  ♦  J  sin  0)  (3) 


3  ■  -rr  (I  cos  0  +  J  sin  0) 


I  sin  0)  .  (4) 


The  radial  and  transverse  components  of  velo¬ 
city  will  be  denoted  by  u  and  v  ,  respectively, 
so  that 

q  -  u(I  cos  0  ♦  J  sin0)+v(J  cos  0  -  I  sin  0  )  ,  (5) 

where  u  «  dr/dt 
and  v  «  r(d0/dt) 

The  differential  equation  governing  the  vari¬ 
ation  with  time  of  the  variables  r,  0,  u  and  v 
is  obtained  by  differentiating  (4)  again  and  sub¬ 
stituting  in  (l).  ffris  results  in  the  equation 


“2  T 

“  r(|^)  (I  cos  0  ♦  J  sin  0) 


cox  9-1  sin  0) 


The  first-order  differential  equation  relating  r 
and  t  is  therefore 

(If 

For  this  differential  equation  to  be  analyti¬ 
cally  integrable,  it  must  be  possible  to  transform 
the  independent  variable  in  such  a  way  as  to  make 
the  right-hand  side  of  (12)  a  perfect  square.  To 
this  end  we  first  express  (12),  in  terms  of  the 
quadratic  roots  of  the  right-hand  side  in  a  dimen¬ 
sionless  form,  as 

(S>!  ■  ■  ■»> 


2  2  2  2 

V  ?  q  -  W  ,  ie  VPQ 
V2(P2  -  Q2)  .  ^  -  U2 


V2-W2  ; 


7  7  f' 

also  P  >  I  >  -  Cf 

and  V2(P2  -  1)0  +  Q2)  -  U2  . 

The  quantity  P  is  always  real  (and  taken  as  posi¬ 
tive)  ,  whereas  Q  has  the  same  character  as  V 
and  is  related  to  the  shape  of  the  orbit:  in  hyper¬ 
bolic  orbits  W  and  Q  are  real;  in  parabolic 
both  are  zero;  and  in  elliptic  orbits  both  are 
imaginary. 

3,  'Minimal'  transformation  of  variable 


whence  the  standard  results 


We  attempt  to  turn  the  right-hand  side  of  (13) 
into  a  perfect  square  by  introducing  a  transforma¬ 
tion  of  the  form 


(S)  •  <£)  ■ « 


It  follows  from  (7)  that 

r  (^r)  “  rv  *  constant  (angular  momentum  per 
unit  mass) , 

so  that 

VR  ,  de  VR 
v  *  —  and 

r  dt  2 


where  a  and  b  are  polynomials  in  n  without 
a  common  factor. 

2  2 

If  Q  ■  -  P  ,  the  right-hand  side  of  (13) 
is  already  a  perfect  square  and  remains  such  under 
any  transformation  of  the  form  (15),  including  one 
found  on  the  assumption  that 

p2  ♦  Q2  i  0  .  (16) 

We  therefore  proceed  with  the  assumption  (16). 
Substitution  of  (15)  into  (13)  yields 
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($S) 


V2(p2b  -  a) (a  ♦  Q2b)/b2  .  (17) 


.  Now,  any  common  factor  of  (a  ♦  Q  b)  and 
(P  b  -  a)  is  also  a  common  factor  of 
[P2(a  +  Q2b)  -  Q2(P2b  -  a) J  and 
[(a  +  Q2b)  +  (P2b  -  a)J  ,  ie  of  (P2  Q2)a  and 
(P2  ♦  Q2)b  .  Thus,  apart  from  a  constant,  there 
is  no  common  factor.  It  follows  that  (a  +  Q2b)  and 
(P2b  -  a)  must  each  be  a  perfect  square.  We  can 
therefore  write 

(a  +  Q2b)  •  [f(n)]2  and  (P2b  -  a)  »  (g (n)  1  2  ,08) 

where  f  and  g  are  polynomials  in  n  without  a 
common  factor. 


In  terms  of  f  and 


2  2  2  2 
P2  ♦  Q2 


.2  .  2 

and  b  »  —z - 

P  +  Q 


09) 


By  substitution  of  (19)  in  (15)  we  obtain  a  trans¬ 
formation  from  r  to  the  new  variable  r\  which 
makes  the  right-hand  side  of  (13)  a  perfect  square, 
whatever  the  form  of  the  functions  f(n)  and 
g(n)  . 

Thus 


pV  -  oV 

f2  ♦  g2 


(20) 


and  then 

m2 

so  that 
dr 


V2(P2  ♦  Q2)2 


.2  2 

J...  L _ _ 

2  2 2 
(f  ♦  g  ) 


_  .  v(P2  .  Q2)  -zfj  • 

Elimination  of  r  between  (20)  and  (21)  yields 


(21) 


dt 

dn 


T  (£g'  -  g£,)  .JA 


<P‘f‘  -  QV>‘  (22) 

where  f'  *  df(n)/dn  g’  "  dg(n)/dn  . 

Then,  without  any  new  assumption  regarding  f  and 
g  ,  (22)  may  be  expanded  in  partial  fractions  as 


-  -2*  (p2  +  q2) (fg’  .  gf 

dn  2W2 


I  1  2u  £g’  -  gf 
2  2  2  2  2 
♦  Qg)  J  VW  P  f  -  Q  g 


(Pf  ♦  Qg) 


(P£  -  Qg) 
2  ' 


(23) 


The  simplest  non-crivial  form  for  f(n)  , 
g(n)  is  linear: 

f  -  F  +  Hn  and  g  -  C  ♦  Kn  ,  (24) 

with  F,  G,  H  and  K  constants.  Then 

fg'  -  gf’  .  FK  -  GH  (constant),  (25) 

and 


The  constants  F,  G,  H  and  K  are  arbitrary;  the 
only  implicit  constraint  on  their  values  is  that 
f  and  g  should  have  no  conmon  factor,  ie  that 
the  expression  in  (25)  should  be  non-zero. 

4.  Kepler's  equation 

With  the  aid  of  (25)  to  (27),  and  using  (14), 
the  equation  (23)  may  be  rewritten  as 

*  -  x 

n/(PH  -  QK)  1  / (PH  ♦  QK)l| 

*L  Pf  -  Qg  Pf  ♦  Qg  Ji 

. (28) 

This  may  be  considered  the  archetype  of  Kepler's 
equation.  When  W  (and  hence  Q)  is  pure  imagin¬ 
ary,  the  logarithm  in  (28)  is  multi-valued,  con¬ 
secutive  values  differing  by  2v i,  ;  thus  identical 
values  of  the  motion  parameters  occur  at  time  inter¬ 
vals  of  2irpt*/W3  .  This  implies  that  the  orbit  is 
periodic,  with  period  T  given  by 

(29) 

Identifying  the  multivalued  logarithm  as  i  times 
an  angle  E  ,  we  can  eliminate  n  in  terms  of  E  . 
Then  we  get 


dt 


■  E  - 


VR 


2(iW) 


?  P(cQ) 


mi 


sin  E 


(30a) 


which  is  a  more  familiar  form  of  Kepler's  equation, 
with  E  the  eccentric  anomaly,  and  with  the  mean 
motion  and  eccentricity  expressed  in  terms  of  V, 

R,  W,  P  and  Q  .  When  W  (and  hence  Q)  is  real, 
on  the  other  hand,  the  logarithm  in  (28)  is  single¬ 
valued  and  can  be  identified  as  the  hyperbolic 
equivalent  of  eccentric  anomaly  (usually  denoted  by 
B)i  elimination  of  n  leads  to  the  hyperbolic 
analogue  of  Kepler's  equation 


VR  P2  +  Q2 
7U2  PQ 


u  , 

7 


(30b) 


The  elliptic  and  hyperbolic  forms  of 
equation  (30)  both  require  the  explicit  determina¬ 
tion  of  the  parameters  P  and  Q  ,  as  defined  in 
(13),  This  is  equivalent  to  finding  the  pericencre 
of  the  orbit,  which  is  a  geometrical  concept 
extraneous  to  the  'minimal1  problem  of  relating  the 
position  and  velocity  at  a  'final'  instant  to  those 
at  an  ’initial’  instant.  To  tackle  the  problem 
without  calculating  P  and  Q  ,  we  can  introduce 
a  constraint  on  F,  G,  H  and  K  ,  so  as  to  set 
n  -  0  at  t  -  0  ,  when  r  -  R  and  u  -  U  .  Then, 
according  to  (20)  and  (21),  we  must  have 

P2F2-Q2G2  -  F2+G2  and  V(P2  +  Q2)FG  -  U(F2  +  G2)  . 

. (31) 

It  follows  that 


4&1  Z.  g,j .  (  (PH  ♦  QK)  /  (Pf  *  Qg) 

p  f  -  q  g 


-(PH  -  QK)/ (Pf  -  Qg) ) / 2PQ 


whilst 


(26) 


P2  -  I 


VFG 

U 


(32) 


d(Pf  ~  S6>  .  (PH  -  QK)  and  d(Pf /■  M  .  (pH  ♦  QK). 
<*1  dn  . (27) 


so  that,  in  terms  of  F  and  G  , 

P2  -  (UG  ♦  VF) /VP  and  Q2  -  (UP  -  VG)/VG  .  (33) 
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to  be  found  from  (42);  then  the  quantities  r,  u, 
v  and  9  can  be  determined  as  follows: 


In  view  of  (14)  and  (33),  F  and  G  are 
constrained  by 

UV(F2  -  G2)  -  (V2  +  W2  -  U2)FG  .  (34) 

The  denominator  in  the  right-hand  sides  of  (20)  and 
(21)  is  f*-  +  g  ;  expressed  in  terms  of  n  , 

f2-g2  -  (F2fG2)  ♦  2(FH  +  CK)n*  (H2*K2)n2  .  (35) 


So  far  H  and  k  are  unconstrained;  it  is 
convenient  to  choose  them  so  that  the  linear  term 
in  (35)  vanishes.  Then  we  put 


H  -  “  ““ 


(36) 


and  get,  with  the  aid  of  (33), 

f2  ♦  g2  -  (F2.G2)  [1  +kV/f21  , 

p2f2  -  Q2g2  -  (F2*G2)  1(1 -KUn/VF)2-  (KWn/VF)2] 

V (P2  +  Q2) fg  -  (F2  +  G2)  [U  +  (V2 +W2-U2)Kh/VF 

-  UK2n2/F2l  . 


. (3?) 

On  substitution  in  (20)  and  (21)  we  get 
R/r  -  [(1  -KUri/VF)2-  (KWn/VF)2]/  [1  .kV/F2! 


-  (U+  (V2  +  W2-U2)Kn/VF-UK2n2/F2]/(!>K2n2/F2]  , 
and  then,  on  denoting  Kn/VF  by  x  » 


and 


Rd  ♦  v2x2) 

Id  -  Ux)2  -  (WX)2] 


(38) 


2  2  2  2  2 
_  dr  .  U  ♦  (V  ♦  W  ~  U^)X  -  UV  X 


1  ♦  v2x2 


(39) 


Then,  on  eliminating  r  between  (38)  and  (39), 
we  obtain 

2  2 

dt  .  2  _ 1  ♦  V  X _ „ 

dx  I)  -  (U  *  WxPTi  -  (U  -  W)xJ  i 

.  “  r  u  ~ w  p  ♦  w  i.  r  x 

^  Ll  -  Ux  *  WX  1  -  Ux  -  WxJ  2w2 


X 


v2  ♦  (U  ♦  W)2 


(U  -  W) 


(1  -  Ux  -  Wx)2  (1  -  Ux  ♦  wx): 


so  that 

dt  -  dji  in ;  -  a  -  ♦  jl 

/ W3  I  -  Ux  ♦  Wx  2W2 

..2 


(40) 


<  f  V2  *  (U  ♦  W)2  .  v2  ♦  (U  -  W)2  ]  ! 

[(U.w)(l  -Ux-Wx'  (U-W)O  -  U  x  ♦  W  x  )| 


and  therefore,  since  t  ■  0  when  x  ■  0  ,  the 
’minimal’  version  of  Kepler's  equation  is 

W3t  _  I  -  Ux  -  Wx  +  RWxfv2*  (U  +  w) 2  .  v2  .  (u  -  w) 2 1 
V  1  -  Ux  ♦  Wx  2u  L  I  -  Ux  -  Wx  1  -  Ux+  Wx  J 

. (42) 

Just  as  for  other  forms  of  Kepler's  equation, 
when  t  is  given,  the  along-track  parameter  x  has 


r  is  given  by  (38)  ; 
u  is  given  by  (39)  ; 
substitution  of  (38)  in  (8)  yields 

v  -  Via  -  Ux)2  -  (Wx)2l/0  +  V2X2)  (43) 

and 

e  -  / (VR/r2)dt  -  /  [2(1  ♦  V2x2)]Vdx 

»  2  arctan  Vx  .  (44) 

Thus  Vx  is  identified  as  the  tangent  of  half  the 
angle  turned  by  the  radius  vector  from  the  initial 
position. 

In  vectorial  form  the  results  can  be  written 
as 

r  -  r  IO  -  V2X2U  +  2VXj)/0  ♦  V2x2) 

-  RtO  -  V2x2)I  ♦  2VXJ)/(0  -  Ux)2*  <Wx)2|  (45) 

and 

q  -  (utO  -  v2x2)I+  2V,'J]  .via  -  v2x2)J- 2VxI)  ) 

*  2  2 
1  *  V  x 

-  (U  +  (W2-U2-V2)x/a  ♦  V2X2))I 

+{V  [1  ♦  (W2  -  U2)x2l/0  +  V2x2))J  .  (46) 

These  formulae  apply  for  all  values  of  U,  V  and 
(non-zero)  R  .  However,  when  W  -  0  (zero  total 
energy,  orbit  parabolic),  equation  (42)  takes  the 
form  0-0  and  has  to  be  expressed  differently, 
eg  by  expanding  both  sides  in  powers  of  W  and 
comparing  the  leading  terms,  allowing  for  the  con¬ 
straint  between  U,  V  and  R  expressed  by  (II). 
The  result  is 

t  .  — ?rx  +  — inL—  .  (47) 

O  -  ux)  30  -  Ux) 

5.  Rectilinear  motion 


Equations  (38)  to  (47)  require  no  modifica¬ 
tion  when  V  ■  0  ,  though  the  right-hand  sides  of 
(43)  and  (44)  become  identically  zero,  expressing 
the  physical  fact  that  the  motion  is  confined  to 
a  straight  line.  Furthermore,  with  the  vanishing 
of  the  terms  containing  V^xZ  in  (39),  that 
equation  becomes  linear  and  allows  the  elimination 
of  v  in  terms  of  the  (radial)  velocity  u  : 


X 


(48) 


Substitution  in  (42)  and  (38)  then  yields  (with 
V  -  0) 


t  - 


and 


(U  +  W)(u  -  W)  RW 
(U  -  W)(u  ♦  W)  "  2u 


(u  -  U) 


X 


u  +  w 
u  -  w 


11 


(49) 


(50) 


7 


In  rectilinear  motion, equation  (II),  defin¬ 
ing  W  ,  reduces  to 


W 


2 


(51) 


and  the  equations  (50)  and  (51)  may  be  combined  as 
r(u2  -  W2)  -  R(U2  -  W2)  -  2u  .  (52) 


Although  this  cannot  occur  in  the  real  world, 
if  the  attracting  central  body  is  reduced  to  a 
point,  then  the  pericentre  of  the  rectilinear 
motion  occurs  at  r  ■  0  .  The  velocity  u  is 
infinite  there  and  changes  sign  from  negative  to 
positive  as  the  point  mass  'bounces  off'  the  centre 
of  attraction.  The  time  from  the  initial  point  to 
the  centre  is  found  from  (49)  by  taking  limits  as 
u  tends  to  infinity: 


U  /  U  +  V  _  RW\  (  . 

lln  im?  —)  ■  (53) 


In  the  further  limit  as  W  -►  0  (rectilinear 
parabola),  we  have 


and 


U  -  u 

V 

RU2  - 

2u  , 

4  /  1 

3  V 

’p) 

4u  . 

2R 

'J? 

"  3U 

(54) 


(55) 


6.  The  Lambert  Problem  and  Lambert's  Theorem 


So  far  we  have  considered  the  'two-point' 
problem  with  the  initial  position  and  velocity  both 
given.  The  results  can  be  applied  also  to  the 
problem  in  which  it  is  the  initial  and  final  posi¬ 
tions  that  are  given,  together  with  the  time  inter¬ 
val,  and  the  velocities  at  the  end  points  are  to 
be  found;  this  is  known  as  the  Lambert  Problem  (or 
the  rendezvous  problem) . 

The  equations  (38)  to  (47)  continue  to  apply, 
but  now  r  and  0  (as  well  as  R)  are  given, 
whilst  U  and  V  (and  hence  W)  are  unknown,  as 
also  is  the  transformed  variable  x  ■  With  t 
given,  equation  (42)  contains  four  unknowns  (x»U,V 
and  W)  constrained  by  (11),  (38)  and  (44). 


The  complexity  of  the  problem  (four  equations 
in  four  unknowns)  is  reduced  if  the  geometry  of  the 
end-point  triangle  (formed  by  the  two  end  points 
and  the  centre  of  attraction)  is  specified  not  by 
the  two  sides  and  included  angle  (R,  r  and  0)  but 
by  the  three  sides,  viz  R,  r  and  c  ,  where  c 
is  the  distance  between  the  end  points,  so  that 


and 


V 


2  2 
X 


2 

c 


2 


-  2Rr  cos  0 


tan^$0 


I  -  cos  0 
I  ♦  cos  0 


C2  -  (R  -  r) 
(R  +  r)2  -  c 


2 

2 


(5ti) 


In  terms  of  R,  r  and  c  ,  the  four  equations  in 
X,  U,  V  and  W  ,  viz  (42),  (38),  (56)  and  (II), 
can  be  written  as 


?! 

,  1  -  «X  -  Wx  .  2WX 

1  -  ♦  w*  (1  -  ux)2  -  (Wx)2 

X  [£w2  .  .  -uxjj  , 

,(57) 

rid 

-  Ux)2  -  (Wx)2]  *  Rd  +  v2x2)  , 

(58) 

[<R  ♦ 

r)2  -  c2]V2x2  -  c2  -  (R  -  r)2 

(59) 

R (U 2  +  V2  -  W2)  -  2u  , 

(60) 

With  the  aid  of  (58)  we  can  eliminate  V 
from  (59)  and  (60),  getting 

t(R  +  r)2  -  c2]  to  -  Ux)2  -  (Wx)2]  -  4R2  (61) 

and 

r!(l -Ux)2- (WX)2]  +  [R(U2-W2)  -  2ulx2  -  R  .(62) 


This  leaves  (57),  (61)  and  (62)  as  three  equations 
in  the  three  unknowns  x»  U  and  W  .  However,  the 
equations  can  be  rewritten  so  that  the  unknowns 
appear  only  in  the  combinations 

W  , 

1(1  -  UX)/R]  ,  «  X  ,  say 
and  [Wx/R]  ,  -  Y  ,  say. 

The  coefficients  in  the  equations  depend  on 
the  known  values  of  u,  t,  c  and  (R  +  r)  ,  but 
are  independent  of  (R  -  r)  ;  thus 


|(R  +  r)2  -  c2]  [X2  -  Y2]  -4  (65) 

and 

(R  +  r)W2[X2  -  Y21  -  2uY2  -  2W2X  .  (66) 

It  follows  that  the  quantities  W,  X  and  Y  in 
(63)  depend  only  on  p,  t,  c  and  R  ♦  r  ,  and  not 
on  R  -  r  .  This  proves  Lambert's  Theorem,  which 
states  that,  with  u  and  energy  (ie  W)  given, 
the  time  of  transit  t  depends  only  on  R  +  r  and 
c  . 


Lambert's  Theorem  permits  the  definition  of 
equivalence  classes  of  triangles  (formed  by  the 
centre  of  attraction  and  the  two  end  points  of  the 
trajectory)  for  which  the  relation  connecting  t, 

W,  X  and  Y  is  identical  for  the  members  of  a 
class;  all  such  members  have  the  same  R  +  r  and 
c  but  are  distinguished  by  the  differing  values  of 
R  -  r  .  This  is  considered  further  by  Gooding6. 

In  particular,  each  equivalence  class 
includes  a  pair  of  degenerate  triangles  for  which 

I R  -  r  |  -  c  , 

ie  such  that  the  two  end  points  and  the  centre  of 
attraction  are  collinear.  The  variation  of  W,  X 
and  Y  with  time  can  easily  be  found  for  the 
degenerate  triangle  (in  which  the  motion  is  recti¬ 
linear),  but  the  results  so  obtained  apply  to  all 
members  of  the  particular  equivalence  class.  We 
distinguish  the  outer  and  inner  end  points  of  the 
equivalent  rectilinear  motion  by  the  suffixes  I  and 
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2,  retaining  the  unsuffixed  notation  for  the  two- 
dimensional  original  problem.  Thus  the  two  sets 
of  radii  and  velocities  are  related  by 


2Rj  *  R  +  r  +  c  ,  2R2  -  R  +  r  -  c  (67) 

and,  in  view  of  (60), 


Ri(“?  '  “2)  *  r2(u2  _  w2)  ■  2»  *  R<y2  +  v2-w2> 


r(u 


2  2 
v  -  V  ) 


.(68) 


From  (49),  the  equation  relating  transit  time 
to  the  other  quantities  is 


t  -  -r  (  In 
w 


<u  ♦  W)(U2  -  W) 

ox.  -  W)(U,  ~  V)*V  ,R2U2-R1U1 


,  (69) 


where  the  skew-symmetric  form  of  the  second  right- 
hand  term  has  been  obtained  with  the  aid  of  (68). 


A  numerical  solution  for  Uj,  U2  and  W  can 
now  be  obtained  using  one  of  the  unknown  veloci¬ 
ties  as  the  iteration  variable  and  calculating 
the  others  with  the  aid  of  (68).  An  accurate 
computing  procedure  for  the  solution,  using  a 
dimensionless  form  of  U1  as  the  iteration  vari¬ 
able,  is  described  by  Gooding6. 


Having  found  the  values  of  U  ,  U2  and  W  for 
the  given  transit  time  t  ,  we  can  determine  the 
other  parameters  in  (63),  since,  from  (48) 

X  -  [1  -U,(U, -U2)/(U2-W2)1/R;  -  [l)u2-w2]/2u 

. (70) 

and 


Y  =  W[(U)  -  U2)/(U2  -  W2)]  /R, 


W(U |  -  U2)/2u  , 
. (71) 


and  therefore,  for  the  individual  members  of  the 
equivalence  class. 


X 

RY 

W 

(72) 

and 

U  - 

w[l 

1  -  RX]  /RY  , 

(73) 

whilst 

from  (36) 

and  (72)  it  follows 

chat 

V  .  |W|  | 

fsL 

-  (R  -  r)2\! 

(74) 

v  tel  1 

u 

*  r)2  -  c2  ) 

Also,  from  the  syrmnetry  of  the  geometry  (or  from 
a  manipulation  of  previously  derived  formulae)  it 
follows  that 


and  hyperbolic  trajectories  (the  parabolic  trajec¬ 
tory  being  a  transitional  case)  and  applying  suit¬ 
able  (and  different)  transformations  of  variable  to 
arrive  at  a  single  algorithm  for  both  types  of 
orbit.  The  transformations  are  introduced  as 
mathematical  tools  of  convenience  and  no  physical 
interpretation  is  offered  by  the  authors  for  the 
end  results.  In  fact,  the  L-B-D  parameters  can  be 
interpreted  in  terms  of  the  Lambert-equivalent 
rectilinear  motion,  as  discussed  in  the  present 
paper.  As  the  L-B-D  parameters  are  all  dimension¬ 
less,  we  first  have  to  express  our  formulae  in  a 
dimensionless  form.  To  this  end  we  must  define  a 
scale  length  and  scale  velocity  related  to  the  geo¬ 
metry  of  the  end-point  triangle  (but  common  to  a 
whole  equivalence  class),  independent  of  the  parti¬ 
cular  time  of  transit  (which  is  the  input  variable 
to  the  problem) .  These  scales  may  be  defined  as 
the  radius  and  the  escape  velocity  at  the  outer 
point  of  the  Lambert-equivalent  rectilinear  triangle, 
ie  _ 


and 


length  scale  =  R} 


whence 


velocity  scale 
time  scale 


A2U/R,) 

/(R?/2«) 


(77) 


The  dimensionless  versions  of  our  formulae 
(68),  (69)  and  (72)  to  (76)  are 

■i 


i  , 

7 


;  In 


•  htt  ■  ■’] 

“  1  , 

(78) 

lu.  +  w]  [  U  -  w]  \ 

'  lu,  -  W][U2  +  Wl  +ulR2U2  Ul>j 

,  (79) 

x  -  R lu ,  -  u2) 

(80) 

-°2«  • 

(8!) 

,  y.„. 

(82) 

R|UI  U2^r(R  +  r) 2 

-c2(  ’ 

[i.w2-u,u2]/[u2 

-v 

(83) 

1  u2  -  (R 

-r>2!! 

(84) 

r'U.  ‘U2l  )(«.  r)2 

-c2(  • 

The  L-B-D  parameters  are  related  to  those  of 
the  present  paper  according  to  the  following  table: 


L-B-D  notation 

Lengths 

s 

c 


Present  notation 


R 


I 


c 


and 

v 

This  ccrople 


u  ■  - 

W[l  ■ 

-  rx]  /rY 

(75) 

-  'Jill 

f  2 

(R  -  r)2V 

(76) 

Irtll 

Ur  + 

tes  the  solution  of  the  Lambert  Problem. 


7,  The  Lancaster-Blanchard-Devaney  method 

4 

Lancaster,  Blanchard  and  Devaney  ("L-B-D') 
developed  a  method  for  solving  the  Lambert  problem 
in  a  coherent  fashion,  by  starting  from  the  clas¬ 
sical  formulae  which  are  different  for  elliptic 


Dimensionless  quantities 

K  (-  q2) 

T 

E 

x 

y  (elliptic) 

(hyperbolic) 


2t 


tW 

W 


z/q 

x-qz 


R2[U2  -  Uj] 


d  (elliptic) 


(hyperbolic) 

f/q  (elliptic) 
(hyperbolic) 

g/q 


(Uj  ♦  W](U2  -  w] 
'!l  ln  Iu|  -  WJ|  U2  ♦  WJ 

(U ,  ♦  w]  [u2  -  w] 

ln  lu,  -  W]  [Uj  +  W| 


tW  [u  j 
w[u, 


U2> 

U2> 


The  multivalued  character  of  the  logarithmic 
term,  when  W  is  imaginary,  is  reflected  in  the 
indeterminate  parameter  m  of  the  L-B-D  analysis. 
No  physical  interpretation  is  offered  for  the 
angular  parameters  a,  B,  y  and  6  ,  which  are 
related2,3  to  the  geometry  of  the  elliptic  or 
hyperbolic  trajectory.  A  final  comment  concerns 
the  separation  of  the  expression  for  the  transit 
time  into  two  terms  described  by  the  L-B-D  function 
$  (and  used  only  for  near-parabolic  trajectories). 
The  two  terms  correspond,  in  the  rectilinear  equi¬ 
valent  motion,  to  the  times  taken  from  either  end 
of  the  trajectory  to  the  instant  of  'bouncing*  off 
the  centre  of  attraction,  the  pericentre  of  the 
rectilinear  motion. 


8.  Conclusions 


The  ’minimalist*  approach  to  orbital  dynamics 
leads  naturally  to  various  basic  results,  in  parti¬ 
cular  as  regards  the  prediction  of  position  and 
velocity  at  a  later  time  when  initial  conditions 
are  given,  and  in  the  solution  of  Lambert's  rendez¬ 
vous  problem.  The  proof  of  Lambert's  Theorem 
becomes  almost  trivial.  Furthermore,  the  elegant 
formulation  of  Lancaster,  Blanchard  and  Devaney 
is  given  a  physical  interpretation. 

This  approach,  applied  to  the  unperturbed 
inverse-square-law  field  of  force,  bypasses  the 
classical  analyses  anchored  to  the  geometrical  pro¬ 
perties  of  the  orbit,  and  avoids  the  complications 
associated  with  the  orientation  of  the  orbital 
conic . 


Orbital  motion  of  a  point  mass  is  described 
by  a  three-dimensional  vector  differential  equation 
of  the  second  order,  and  therefore  requires  six 
parameters  (orbital  elements)  to  describe  it  fully. 
Five  of  these  may  be  used  to  describe  the  path  in 
space,  and  the  sixth  to  relate  it  to  absolute  time 
(as  by  stating  the  time  at  pericentre  or  the  cros¬ 
sing  of  a  reference  plane).  The  classical  set  of 
spatial  elements,-3  consists  of  two  describing  the 
shape  and  size  of  the  orbit  (a,e)  and  three  angles 
giving  the  orientation  of  the  axes  of  the  orbit 
(i,ft,ui).  An  alternative  spatial  set  is  provided  by 
the  five  independent  components  of  the  (mutually 
orthogonal)  angular  momentum  and  eccentricity 
vectors,  which  can  be  defined  in  terms  of  ine  posi¬ 
tion  and  velocity  vectors.  In  principle  it  does 
not  matter  which  six  elements  are  used  to  describe 
the  orbit,  provided  they  are  not  subject  to  an 
internal  constraint.  In  practice  some  choices 
make  the  problem  ill-conditioned  (eg  pericentre- 
related  parameters  in  near-circular  orbits,  or 
node-related  parameters  in  near-equatorial  orbits). 
Additional  difficulties  may  arise  in  numerical 
work,  with  loss  of  accuracy  when  small  differences 


between  large  quantities  have  to  be  calculated. 

There  can  therefore  be  no  hard-and-fast  rule  which 
would  state  that  a  particular  set  of  orbital  elements 
is  better  than  all  others  in  all  circumstances.  The 
success  achieved  here  with  the  direct  relationships 
between  the  position  and  velocity  vectors  at  two 
particular  points,  supplemented  only  by  the  velocity 
W  ,  suggests  that  perturbation  effects  might  be 
expressible  in  terms  of  the  (slow)  variation  of  the 
position  and  velocity  at  a  selected  epoch. 
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introduced  by  other  authors  through  a  purely  mathematical  analysis.  6r  >c  .  •  t.  )£— 
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